Larger Domains from Resonant Decay of Disoriented Chiral Condensates 
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The decay of disoriented chiral condensates into soft pions is considered within the context of a 
linear sigma model. Unlike earlier analytic studies, which focused on the production of pions as 
the sigma field rolled down toward its new equilibrium value, here we focus on the amplification of 
long-wavelength pion modes due to parametric resonance as the sigma field oscillates around the 
minimum of its potential. This process can create larger domains of pion fluctuations than the usual 
spinodal decomposition process, and hence may provide a viable experimental signature for chiral 
symmetry breaking in relativistic heavy ion collisions; it may also better explain physically the large 
growth of domains found in several numerical simulations. 



PACS 25. 75. +5 Preprint HUTP-97/A107, iep-ph/9801307 



Experiments at the Relativistic Heavy Ion Collider at 
Brookhaven and at the Large Hadron Collider at CERN 
may soon be able to probe many questions in strong- 
interaction physics which have until now been studied 
only on paper or simulated on a lattice. One major area 
of study concerns the QCD chiral phase transition. In 
relativistic heavy ion collisions, it is possible that non- 
equilibrium dynamics could produce "disoriented chiral 
condensates" (DCCs), domains in which a particular di- 
rection of the pion field develops a non-zero expecta- 
tion value. |l]-[|] These domains would then decay to the 
usual QCD vacuum by radiating soft pions. Preliminary 
searches for DCCs by the MiniMax Collaboration in pp 
collisions at Fermilab have thus far not found evidence 
for the production and decay of DCCs Q , though they 
are far more likely to be created in upcoming heavy ion 
collisions. Thus, understanding their possible formation 
and likely decay signatures in anticipation of further ex- 
perimental work is of key importance. 

If these domains grow to sufficient size (on the order of 
3-7 fm), such an experimental event would be marked by 
a particular clustering pattern: some regions within the 
detector would measure a large number of charged pions 
but few neutral pions, while other regions of the detec- 
tor would measure predominantly neutral pions with few 
charged pions. || Defining 1Z to be the ratio of neutral pi- 
ons to total pions, 1Z = n OT o/(n 7r o+n 7r +-r-n 3r -), it has been 
demonstrated that the probability for measuring various 
ratios K in DCC events obeys P(ll) = (4ft)- 1 / 2 , which, 
especially for small- 1Z, may be easily distinguished from 
the isospin-invariant result of P(1Z) — > 6(1Z — 1/3). ||] 
(Detecting the decay of such DCCs could be improved 
by measuring the two-pion correlation functions Q , and 
from enhanced dilepton and photoproduction Q , in addi- 
tion to studying the fraction of neutral pions produced.) 
The production and subsequent relaxation of such DCCs 
may also explain the so-called "Centauro" high-energy 
cosmic ray events, in which very large numbers of charged 
pions are detected with only very few neutral pions. [g|J^] 

However, as emphasized in if the disoriented do- 
mains do not grow to such large scales within heavy ion 
collisions, such experimental signatures become less and 



less easy to distinguish from the isospin-invariant case. 
Even if DCCs are produced following a heavy ion colli- 
sion, if the domains do not grow to be "large" (that is, 
several fm), then the detector would sample so many of 
these discrete domains within a given run, each of which 
with the pion field aligned along some random direction, 
that the clustering effects would be washed out. A cru- 
cial question, then, is whether or not sufficiently large 
domains might grow in the noncquilibrium aftermath of 
a heavy ion collision. 

Several authors have considered the amplification of 
long-wavelength pion modes from the decay of DCCs in 
the context of a linear sigma model. [^]j|,||,|9| The rele- 
vant degrees of freedom are modeled by the scalar fields 
a and 7? , which may be grouped together as $ = (cr, 7?) . 
Above the critical temperature, when the system is in the 
chirally-symmctric state, the effective potential for $ is 
0(4) symmetric, and ($) = 0. To model a strongly- 
nonequilibrium situation, Rajagopal and Wilczek con- 
sidered a quench: as the quark-gluon plasma in the 
interaction region between the colliding nuclei expands 
and cools, the effective temperature may fall quickly 
to T <C T c . Because the zero-temperature potential 
is not chirally-symmetric, domains form, and it takes 
some time for the fields to evolve from (cr) = (7?) = 
to the new equilibrium values, (cr) ^ 0, (7?) = 0. Fol- 
lowing the quench, the fields relax to these new equilib- 
rium values according to the effective potential, V($) = 
j(& 2 — v 2 ) 2 — Ha; that is, the a field 'rolls down' from 
a ~ to a ~ v. 

Numerical simulations [^]|| reveal a large amplification 
of soft pion modes from the relaxation of the nonequi- 
librium plasma. Previous authors have attempted to 
explain these numerical results analytically in terms of 
spinodal decomposition: during the time that a rolls 
down toward v, pion modes with wavelengths satisfying 
k 2 < X(v 2 — ( < I >2 )) will grow exponentially. However, 
under the usual quench scenario, the time it takes for 
a to roll to v, and hence the maximum domain size for 
the DCCs, remains too small to produce clear experi- 
mental signatures. Under this scenario, domains typi- 
cally remain pion-sized, ~ 1.4 fm. (See, e.g., This 
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physical mechanism alone therefore remains incapable of 
explaining the large domains found in numerical simula- 
tions. 

Building on earlier work in |l(i[| , we consider here a 
physically distinct process which could produce larger 
domains of DCC, and hence might better explain the 
significant clustering observed in numerical simulations. 
Rather than amplification of pion modes while a rolls 
down its potential hill, we focus on the parametric am- 
plification of pion modes as a oscillates around the min- 
imum of its potential. Because this is a distinct process, 
the growth of domains due to parametric resonance, un- 
like the growth of domains due to spinodal decomposi- 
tion, may reach scales on the order of 3-5 fm. 

This means of DCC decay is similar to cosmological 
post-inflation reheating. An early attempt was made 
to apply the reheating formalism of jy] to the decay 
of DCCs due to parametric resonance in |!o| . However, 
the analytic tools for studying the noncquilibrium, non- 
perturbative dynamics of such resonant decays have im- 
proved since this early work on reheating, and the ear- 
lier approximations, while at times qualitatively informa- 
tive, prove quantitatively unreliable. Most important, 
this earlier study jHj] approximated er's oscillations as 
purely periodic, in which case the equation of motion for 
the pionic fluctuations reduces to the well-known Math- 
ieu equation. Ignoring the nonlinear, anharmonic terms 
(such as Act 4 ) in the evolution of a then yields the pre- 
diction of an infinite hierarchy of resonance bands, with 
decreasing characteristic exponents. Yet given the non- 
linear equation for cr, the equation of motion for the pi- 
ons reduces instead to a Lame equation, which, in the 
cases of interest, has only one single resonance band, with 
a different value for the amplified modes' characteristic 
exponent. As emphasized in |12| ], these two differences 
combined can change dramatically the predicted spec- 
tra from parametric resonance; to be useful in making 
contact with experiments, these nonlinearities must be 
attended to, as in the present study. (Furthermore, the 
authors of |[o| did not consider the size of domains cre- 
ated by the parametric resonance, as considered here.) 
Instead, we draw on the more recent studies of reheating 
in [p~2| p~5]| to consider the question of DCCs and their 
resonant decay. 

Following we consider a quench scenario: the tem- 
perature of the plasma drops quickly from above the crit- 
ical temperature (with ($) = 0) to near zero. The effec- 
tive Lagrangian density following the quench is given by 

£ = -\{d^) 2 ~ ± (&-v 2 ) 2 + H*. (1) 

Here H is an external field which breaks the chiral sym- 
metry and picks out the a direction as the true minimum. 
The pion mass is proportional to H . The true vacuum is 
characterized by ($) = {f n ,0), where = 92.5 MeV is 
the pion decay constant. In the limit as H — > 0, — > v . 
In the following, we neglect H in the resulting equations 
of motion, but add by hand a pion mass — 135 MeV; 



we also set A = 20.0 and v = 87.4 MeV, which yield 
TO(T = (2A/2 + ml) 1 / 2 = 600 MeV. These standard val- 
ues for the parameters are chosen, as in to fit 
low-energy pion dynamics. 

As a first approximation, we neglect effects due to the 
expansion of the plasma. Obviously the expansion of the 
plasma plays a crucial role, at least for early times fol- 
lowing the collision, in dropping the temperature below 
the critical temperature. (Some work has been done to 
incorporate analytically the effects of cosmological ex- 
pansion in the resonant decay of a massive inflaton [jl4| , 
which may be useful in improving future analytic stud- 
ies of DCCs and their decay). We also ignore noise and 
other medium-related effects on the resonance; as demon- 
strated in the context of post-inflation reheating, such ef- 
fects do not generally destroy the parametric resonance, 
but rather enhance it. |]l6| 

We study the nonequilibrium, nonperturbative dynam- 
ics by means of a Hartree approximation, by writing 
c(t,x) = <Jo(t) + 6a(t,x), and replacing if 3 — > 3(tt 2 } tt 
and 7? 2 — ► (7F 2 ). The vacuum expectation value may 
be written in terms of the field's associated (Fourier- 
transformed) mode functions as (if 2 ) = J d 3 k\irk\ 2 / (2ir) 3 . 
Because the 5a fluctuations decouple from the pion 
modes in this approximation, we will focus below on the 
pionic fluctuations. 

Within a given DCC domain, the pion field will be 
aligned along some particular direction, n ff , in isospin 
space. We will therefore writ e if — \ n n . In ter ms o f the 
dimensionless variables r = V Xv 2 t and n = k/ V Aw 2 , and 
the scaled field tp(r) = a (t)/v, the coupled equations of 
motion take the form: 

<f" + (<y5 2 - 1 + £,) ¥> - 0, 

Xk + (P 2 + f 2 + V*) Xk s 0, (2) 
where primes denote d/dr, and we have defined 

M = m^/\/Aw 2 , p 2 = k 2 + M 2 - 1, E T = (n 2 )/v 2 . (3) 

Note t hat with the values of the parameters assumed 
here, VXv 2 = 390.9 MeV, and M 2 = 0.12. These equa- 
tions of motion are conformally equivalent to those for 
massless fields in an expanding, spatially-open universe, 
and hence we may apply the techniques of ]l5{] to study 
their solutions. 

We are interested in the growth of 7? modes as <tq oscil- 
lates around v. Having begun, following the quench, near 
<7o ~ 0, (To will roll down its potential hill toward v. The 
rolling field will at first overshoot the minimum at v, and 
then begin oscillating around v. The amplitude of these 
oscillations will eventually be damped by the transfer of 
energy from this oscillating zero mode into the 7r fluctua- 
tions. For early times after these oscillations have begun, 
however, the amplitude of cto will remain nearly constant. 
In this strongly-coupled system, unlike in the weakly- 
coupled inflationary case, the a field will execute only a 
few oscillations before settling in to its minimum. Yet, 
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as we see below, even these few oscillations could prove 
significant, since most particle production via paramet- 
ric resonance occurs in highly non-adiabatic bursts, when 
the velocity of the oscillating field passes through zero, 
jnj Furthermore, because the system has been quenched 
from its initial, chirally-symmetric state, we assume that 
Sjr is small at the beginning of ao's oscillations. (The 
fact that spinodal decomposition alone cannot produce 
large DCC domains is equivalent to remaining small 
while a rolls down its potential hill.) Then we may solve 
the coupled equations for early times after the oscilla- 
tions have begun, and study the growth of the fluctua- 
tions 7?fe . Because Co begins oscillating quasi-periodically, 
certain pion modes will be amplified due to parametric 
resonance. 

The resonance will fade once the backreaction term, 
S,r, grows to be of the same order as the tree- level terms, 
such as ip 2 . To study the behavior of the pionic fluctu- 
ations, we solve the coupled equations of (0) for early 
times after the beginning of cto's oscillations, when E w 
may be neglected. Thi s lasts up to the time r en( j, deter- 
mined by S,r(r en d) = p 2 (r), where an overline denotes 
time-averaging over a period of ^'s oscillations. 

Assuming that a^s oscillations begin once <7o reaches 
its inflection point, a- mn = it will roll past the min- 

imum and up to the point at which V{a) = T^(<7i n fl), be- 
fore rolling back down through v. This sets ip = ^/5/3. 
Because this definition of the initial amplitude is some- 
what arbitrary, we study the resonance effects for 1^0 in 
the range 1 <J£o < V%- In the range f < tp < V2, ¥>(t) 
oscillates as Ifl2l 



V{ T ) = <Po dn (jt, v) 



(4) 



where dn(it, v) is the third Jacobian elliptic function, 7 = 

(fo/VZ, and v = ^2(1 - tp^ 2 ). Eq. (@) holds for T < 
Tend- The dn- function oscillates between a maximum at f 
and a minimum at (f — v 2 ) 1 / 2 , with a period of 2K{v)/^), 
where K(v) is the complete elliptic integral of the first 
kind. @ 

With ip(r) oscillating as in Eq. (Q), the equation of 
motion for \k becomes the Lame equation of order one. 
A solution for the pion modes Xk (t) may thus be written 
in the form fflM: 



(5) 



C/f ) (r)=A(r)e X p (±/i K (^) 7 r) 



Here A(t) is a periodic function, normalized to have 
unit amplitude, and n K {y) is the characteristic expo- 
nent (also known as the Floquet index). The form 
of p K depends on both k and v. Clearly, whenever 
Re [/ji k (v)} ^ 0, the coupled modes will be exponentially 
amplified. The exact relation between the Uk modes and 
Xk (and hence 7?) depends on the assumed initial condi- 
tions following the noncquilibrium quench. If we make 
the usual assumption, that Xfc( T o) = 1/ y / 2w p (to) and 
X ' k (r ) = -i^/lu p (t )/2, with to 2 (t) = p 2 + ^ 2 (r) ffJQ, 



then the Xk modes may be written as a linear combina- 
tion of U ( k ±} as in ||^||. 

In |l5|], these coupled equations were studied for the 
range tpo > \/2. Proceeding in exactly the same way, 
solutions may be found for the case 1 < <po < 
The characteristic exponent has non-zero real parts only 
within a single resonance band, given by 



(i _ AT 2 ) - yi- p 2 {2-p 2 )<n 2 
-M 2 ) + yi-p 2 (2-^). 



< 



(6) 



The resonance band includes modes with k < for all 
values of </?o > 1.23, that is, even fo r am plitudes of the os- 
cillating field smaller than tpo — ^/5/3. As in jTJ[, n K {v) 
may be written in terms of a complete elliptic integral of 
the third kind. The real part of Uniy) is plotted in Fig. 1. 
Near the center of the resonance band for a given value 
of Re [/i K ] ~ 0.1 — 0.3. Note that as in the numerical 
simulations of ^|J|], the strongest amplification (indicat- 
ing greatest particle production) occurs for k < m^. The 
maximum values of Re [fi K ] fall in the k — ► limit. 




FIG. 1. Re[/j, R (v)] as a function of both the dimensionless 
momentum, k, and the dimensionless initial amplitude of ao's 
oscillations, ifo. In these units, = 0.35; the largest expo- 
nents, and hence the strongest resonance, occur for k < m^. 

Given Re\p, K ], one can determine r cn d, based on the 
growth of S^. If T on( j is large enough, then observable 
domains of DCC could be formed and detected. Within a 
given domain, E^r) is given as an integral over |xfe(r)| 2 . 
To evaluate r on d, we solve numerically the equation 



(Tend) 



res. band 



2tt 2 



|Xfc( r ond)| 2 = (fi 2 , 



(7) 



where the integral extends over the single resonance 
band. The time-average of the (dimensionless) oscillat- 
ing cr field, ip 2 , may be written in terms of -E(^), the 
complete elliptic integral of the second kind, using the 
integral of dn 2 (u, v) over a period of its oscillations (see 
]l7|). Eq. ([?]) then yields t C nd ~ 5 fm/c over most of 
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the range 1 < ipo < \/2. This should be compared with 
the usual spinodal decomposition scenario, in which the 
pionic fluctuations would be exponentially amplified only 
for the brief period i S pinodai — V2/m a ~ 0.47 fm/c. 

In order to find the characteristic sizes to which DCC 
domains may grow before the parametric resonance is 
damped, we may follow |2Cj and evaluate the two-point 
correlation function: 

r ffih. 

D(t, r) = (n(t, r)n(t, 0)> = J —^\n k f 

= (^)- 3/2 f ^ioM| Xfe (r)| 2 , (8) 

with the dimensionless length defined by x = %/ \v 2 v. We 
have also used 7i>(t) = Xk(j) within a given domain. 
We may solve this integral in the saddle-point approx- 
imation, making use of Eq. (11.4.29) of |17| , with the 
result that 

D(t, rJocE^tJexpf-rV&W). 
e D (t) = 4 7 |i?e {d 2 ^( V )/dn 2 ) | Kmax (i/VA^). (9) 

Eq. (Q) reveals that the domain size grows as t 1 / 2 . The 
maximum correlation length, £,o(icnd), is plotted in Fig. 
2. 
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FIG. 2. Maximum domain size for DCCs, ^D(tend), as a 
function of ifo, in units of fm. 

Over much of the range 1 < ipo < \/2, ^c(iend) lies 
between 3-5 fm. The usual spinodal decomposition pro- 
cess, in the absence of the 'annealing' studied in ||, on 
the other hand, can only create domains of order 1.4 fm. 
||. If 'annealing' is effective, domains from spinodal de- 
composition may grow to 3 — 4 fm; yet independent of the 
dynamics as er rolls down its potential hill, it still must 
end its evolution by oscillating around v, and the reso- 
nant production of pions would follow as studied here. 

Parametric resonance offers a promising means of pro- 
ducing observable signals from the production and decay 
of disoriented chiral condensates in the aftermath of rel- 
ativistic heavy ion collisions. Unlike spinodal decompo- 
sition alone, this physical process may explain the sign- 
ficant growth of DCCs found in numerical simulations. 
Future analytic studies should better include effects from 
the expansion of the plasma, and from 8<r — tt and 7? — 7? 
scatterings, which are neglected here in the Hartree ap- 
proximation. If these effects remain subdominant, how- 
ever, then the resonant decay of DCCs should produce 



low-momentum pions with a distribution observably dis- 
tinct from the isospin-invariant case. 
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